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ABSTRACT 

CN '. 

£SJ , If primordial fluctuation is Gaussian distributed, higher-order cumulants of the 

cosmic fields reflect nonlinear mode coupling and provide useful information of 
00 ! gravitational instability picture of structure formation. We show that their expected 

deviation (sample variance) from the universal values is nonvanishing even in linear 
' theory in the case where observed volume is finite. As a result, we find that the relative 

sample variance of the skewness of the smoothed velocity divergence field remains as 
2 ■ large as ~ 30% even if the survey depth is as deep as ~ 150/i _1 Mpc. 

Subject headings: cosmology: theory — large-scale structure of the universe 
«J ■ 1. Introduction 

d; 

It is now commonly believed that the large-scale structure in the Universe observed today 
is explained by gravitational evolution of small initial density inhomogeneities (Peebles 1980). 
These initial fluctuations are usually assumed to obey random Gaussian distribution which is 
not only plausible from the central limit theorem but is also predicted by standard inflation 
models (Guth &; Pi 1982, Hawking 1982, Starobinski 1982). Quantitative analysis of statistical 
measures of present-day cosmic fields are very important to confirm or disprove the structure 
formation scenario based on gravitational instability from primordially Gaussian fluctuations. 
Their higher-order cumulants provide useful tools for this purpose. 

In addition to the anisotropy of the cosmic microwave background radiation (CMB), 
distribution of galaxies and the peculiar velocity field are basic measures for the statistical analysis 
of the large-scale inhomogeneities in the Universe. The former has been widely investigated 
and there are two ongoing large redshift surveys now, namely, the Anglo-Australian 2dF Survey 
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(Colless 1998) and the Sloan Digital Sky Survey (Gunn & Weinberg 1995) which are expected to 
revolutionarily improve our knowledge of three-dimensional galaxy distribution. 

We should notice, however, that what we can directly observe is the distribution of galaxies 
whereas what we can discuss from the first principle is the distribution of the underlying matter. 
In spite of the rapid increase of observational data, our understanding of the relation between 
distribution of galaxies and that of underlying gravitating matter, namely biasing (Kaiser 
1984), is far from satisfactory. This hampers straight forward comparison between theories and 
observations. 

In contrast to the number-density field of galaxies, the cosmic velocity field reflects the 
dynamical nature of underlying matter fluctuation and is basically independent of the poorly 
understood biasing relation at least on large scales (Dekel 1994). This is a fundamental merit of 
the cosmic velocity field. On the other hand, we must point out that the survey depth of the comic 
velocity field is currently limited only to L ~ 70/i _1 Mpc around us even in the case of recent Mark 
III Catalog of Galaxy Peculiar Velocities (Willick et al. 1997), which is much smaller than that of 
the redshift surveys^]. Therefore uncertainties due to the finiteness of our survey volume, or the 
sample variance, must inevitably become large and are therefore very important in the analysis of 
the cosmic velocity field. 

The higher-order cumulants of velocity divergence field have been extensively investigated in 
the framework of nonlinear perturbation theory, and are expected to work as useful quantities 
in observational cosmology, for example, to constrain the density parameter independent of the 
biasing (Bernardeau et al. 1995, 1997). In this Letter we investigate the sample variance of the 
higher-order cumulants of the velocity divergence field assuming that the initial fluctuation obeys 
isotropic random Gaussian distribution. Our formalism is also applicable to the density field and is 
similar to Srednicki (1993) who analyzed the skewness parameter of cosmic microwave background 
radiation. 

In the previous Letter (Seto & Yokoyama 1998) we discussed the sample variance of the 
second-order moment or the variance of a component of the peculiar velocity and that of the 
linear density fluctuation, whose expectation values and sample variances are of the same order 
in perturbation. In the case of higher-order cumulants, their expectation values vanish in linear 
theory and are generated from nonlinear mode coupling in higher-order perturbation theory. 
Nevertheless their sample variance is nonvanishing even at the linear order. Thus the sample 
variance is expected to be much more important for higher-order cumulants. In this Letter we 
compare the expectation values of the lowest-order contribution of the higher-order cumulants of 
these fields and their sample variances predicted by linear theory. 

It is true that there are other sources of errors in the observational analysis of these fields. 



1 Bernardeau et al. (1995) used L ~ 40h 1 Mpc as a practical current limit of high-quality data of the peculiar 
velocity field. 
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Using Monte Carlo calculations we could basically take various effects into account at one time 
(e.g. Borgani et al. 1997). But the sample variance due to the finiteness of the survey region can 
be regarded as a fundamental limitation in the sense that this uncertainty is independent of how 
accurately we could measure the cosmic fields in a specific region in the Universe. In addition, 
to investigate the sample variance by means of a numerical simulation, we need a simulation box 
much larger than the (mock) survey region, as the sample variance is heavily weighted to Fourier 
modes which are comparable or greater than the survey depth. Considering these two factors, 
the simple analysis presented in this Letter is a useful and convenient approach to estimate a 
fundamental limitation in the observational determination of higher-order cumulants of cosmic 
fields. 



2. Formulation 

We denote the density contrast field by <5(ai) and the velocity divergence field by 
9{x) = H^V ■ V(x) where V(x) is the peculiar velocity field and Hq is the Hubble parameter. 
At the linear order, which is indicated by the suffix "lin" hereafter, we have the following relation. 

Mas) = -/MM*). (!) 

where the function / is the logarithmic time derivative of the logarithm of the linear growth rate of 
the density contrast 5n n (x) and is well fitted by /(fio) — ^o' 6 with Qq being the density parameter 
(Peebles 1980). We define the linear dispersion of these fields as 

a 2 ee (Sl n (x)) , ^(te), (2) 

where (X) represents to take an ensemble average of the a field X. From equation (1) the linear 
root-mean-square (RMS) fluctuation of the velocity divergence field a$ is written in terms of a 
and f(flo) as 

a e = f(n )a. (3) 

These two quantities a and a$ work as the expansion parameters for perturbative analysis in this 
Letter. 

In the observational study of the cosmic fields in the framework of perturbation theory a 
smoothing operation is crucially important to get rid of strong nonlinearities on small scales and 
noises due to the discreteness of galaxies which work as tracers of these fields. In this Letter we 
only discuss fields smoothed with a Gaussian filter defined by 

Here R s is the smoothing radius but we omit its explicit dependence in most part of this Letter 
for notational simplicities. 
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Next we briefly summarize the expectation values of the third- and forth-order cumulants 
for two fields, 5(x) and 9(x). We introduce their first-nonvanishing contributions predicted by 
higher-order (nonlinear) Eulerian perturbation theory (e.g. Peebles 1980). First we assume that 
the power spectrum of density fluctuation has a power-law form characterized by a single power 
index n (> —3) as 

P(k) cx k n . (5) 

In this case the third-order cumulants or the skewness of S(x) and 9(x) smoothed with the 
Gaussian filter @ are evaluated perturbatively as follows. 

5(xf) = S 3 (n)a* + 0(a 6 ), (9(xf) = S 36 (n)4 + 0(a 6 e ), (6) 

where the factors S3 (n) and 630(71) are of order unity and have been given by Matsubara (1994) 
and Lokas et al. (1995) in terms of the hypergeometric function F. 

_ + 3 n + 3 3 1\ / 8\ /n + 3 n + 3 5 1\ 

S3{n) = 3F I— — 2' I)" [ n+ 7 r l— — 2'4 ' (7) 



639(71) = 



mo) 



3p n + 3 n + 3 3 IN / 16N /n + 3 n + 3 5 1 



2' 4/ V 7 7 V 2 ' 2 ' 2' 4 



(8) 



Here we neglect extremely weak dependence on cosmological parameters except for the function 
/(fio) in Sse(n). In principle, equations (0) and (§) are valid only for a pure power-law spectrum 
as equation (||) but we extrapolate them to more realistic power spectra with an effective power 
index n(R s ) defined at the smoothing scale by the following equation (see Bernardeau et al. 1995): 

n{Rs) ~ " 3 " dl nj R s " " 3 " dlni? s • (9) 

In the same manner the forth-order cumulants, or the kurtosis, of S(x) and 9(x) are written 
perturbatively as follows. 

'<5(a0 4 - 3cr 4 ) = S 4 (n)a 6 + 0(a s ), (10) 
9(x) i -2,af) = S 4 e(n)a$ + 0(a$). (11) 



Lokas et al. (1995) derived analytic formulas for 64 (n) and S40 (n) based on higher-order 
perturbation theory and evaluated them numerically for —3 < n < 1. In the present analysis we 
use fitting formulas for 64(77) and 649(77) given in their paper. 

In observational cosmology, we usually estimate an ensemble average (X(x)) of a field X by 
taking its volume average A(X, V) in an observed patch V, 

A(X, V) = i J X(x)d 3 x => (X) . (12) 

We can commute the ensemble average with a volume integral above to obtain 

(A(X, V)) = (X) . (13) 
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Thus the ensemble average of the volume average A(X,V) is identical to the universal value (X). 
However, the observed value A(X, V) in one specific patch V is expected to fluctuate around its 
mean (X) because of the spatial correlation and inhomogeneity of the field X(x) beyond the 
patch V. These fluctuations are nonvanishing even in linear theory and we define its RMS value 
E\i n (X, V) as follows. 

E liQ (X, V) = ({A(X nn , V) - (X lin )} 2 ) 1/2 . (14) 

Our basic strategy is to compare the magnitude of this linear sample variance E^ n (X, V) with 
the expectation value {X). This fluctuation should be smaller than the expectation value (X); 
otherwise the particular value of A(X, V) obtained in one survey volume would lose its universality 
and one could not extract any cosmological information from it. 

Let us now calculate the linear fluctuation En n (X, V) for the skewness and the kurtosis of the 
velocity divergence field. Using the nature of the multivariate Gaussian variables, we obtain the 
sample variance of the skewness of 9{x) as 

EL(0 3 ,V) = M J^x J v d 3 yZ(r xy ){3 + 2Z(r xy ) 2 }, (15) 

where we have denoted the separation between two points x and y by r xy = \x — y\ and defined 
the normalized linear two-point correlation function H(r) as 

w , (0Ux)e Vm (y)) _ (5 Vm (x)6 Vm (y)} _ [°° k 2 dk sin kr 2 2 

Ar*y)- - ~ 2 -J o 27r2(J 2 kr nfc)exp(-fc K s ). (lb) 

In the same manner the linear fluctuation for the forth-order and second-order cumulants are 
given as follows. 

EKe'-SatV) = -^j v d 3 xj v d 3 yZ(r xy ) 2 {3+~(r xy ) 2 }, (17) 
EL(0 2 ,V) = ^J^xJ v d 3 y~(r xy ) 2 . (18) 

In the next section we calculate the ratio of E\[ n (X, V) to the lowest nonvanishing order of 
(X). Here one should notice that the skewness and kurtosis given in equations (6), (10), and (11) 
are obtained from higher-order contributions in perturbation in contrast to En n (X, V) obtained in 
linear theory. We summarize the order of the expansion parameter gq for these ratios below. 

E(d 2 ,V) 



(0 2 ) 
E(6 3 ,V) 

(0 3 ) 



0(1), (19) 



0(a e % (20) 



EjO'-^lV) _ ,- 2) (21) 



So far we have mainly discussed linear fluctuation of the velocity divergence field 9{x) but 
extension to the case of the density field 5(x) is simple and straightforward (see eq.|l6|). 
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3. Results 

In this section we calculate the ratio E\ in (X,V)/\ (X) | using specific cosmological models. For 
calculational simplicity we assume that our survey patch V is a sphere with radius L and volume 
V = (47r/3)L 3 . In this case we can simplify the six-dimensional integral of equations (|i~5|), (|i~7l), 
and (|j~8| ) to a three-dimensional one owing to the rotational symmetry (Seto & Yokoyama 1998). 

We investigate two cold-dark-matter (CDM) models with different density parameter Qq, 
namely, Slo = 0.3 (open model) and = 1-0 (Einstein de-Sitter model) both with vanishing 
cosmological constant^ and the Hubble parameter h = i7o/(100km/sec/Mpc) = 0.7. As for the 
initial matter fluctuation we use CDM power spectrum given in Efstathiou et al. (1992) as 

, _ . -B k , . 

P ( k ) = 7 ( 22 ) 

where a = (6.4/r)/ i - 1 Mpc, (3 = (3.0/r)/i~ 1 Mpc, 7 = (1.7/r)/i- 1 Mpc, fi = 1.13, and the 
normalization factor B = (96tt 2 / 5)^Iq 1 ' 54 Hq 4 (Q rms /Tq) 2 with the current temperature of CMB 
To = 2.73K and the quadruple fluctuation amplitude of it Q Tms = 15.3/iK from 4yr COBE data 
(Gorski et al. 1996). We fix the shape parameter r by T = h£lo. 

As explained before, the cosmic velocity field is considered to be less contaminated by the 
poorly understood biasing effect but its survey depth is much smaller than that of the redshift 
surveys of galaxies. Therefore we mainly consider a typical observational situation of the cosmic 
velocity field and adopt a Gaussian filter with R s = 12/i _1 Mpc following the POTENT analysis 
(Bertschinger & Dekel 1989, Dekel 1994). 

Using the formulas given in the previous section, we plot the sample variance due to the 
smallness of the survey volume in Fig.l. The expansion parameter (a, ag) is (0.37,0.18) for the open 
model flo = 0.3 and (0.38,0.38) for the Einstein de-Sitter model f^o = 1-0 in our case. We have 
0{x) oc 5(x) at the linear order and thus the relative fluctuations of the second-order moments are 
identical for these two fields (thick solid lines in Fig.l). For the higher-order cumulants, nonlinear 
mode coupling arises in a different manner for 9(x) and 5{x) and their relative fluctuations are no 
longer identical. 

As is seen in Fig.l, at the current survey depth L ~ 40/i -1 Mpc (Bernardeau et al. 1995), the 
sample variance of the skewness of the velocity divergence field (#(a:) 3 ) remains as much as the 
expectation value itself. To reduce it smaller than ~ 30%, we have to take the survey radius L as 
deep as L ~ 180/i _1 Mpc for the open model and L ~ 150/i _1 Mpc for the Einstein de-Sitter model. 
These values are much larger than the current observational limit. Bernardeau (1995) proposed 
a method to estimate the density parameter using a relation (0 3 ) / (6 2 ) cx Qq ' 6 . Analysis above 



2 Cosmological constant is almost irrelevant in our analysis, even if it is nonvanishing within the current 
observational limit. 
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show that even if we could take the survey radius L as big as ~ 150/i Mpc, the error of such 
estimation of Oo would remain as large as ~ 50%. 

4. Summary 

In this Letter we have discussed the magnitude of the sample variance in the observational 
determination of the reduced higher-order cumulants of smoothed density and velocity divergence 
fields. We have compared the sample variance predicted by linear theory with the lowest-order 
nonvanishing contribution to the cumulants assuming that the primordial fluctuation was random 
Gaussian distributed. We have paid much attention to the velocity divergence field as (i) it is less 
contaminated by the biasing relation and extensively investigated in the framework of perturbation 
theory but (ii) velocity survey is currently limited to a relatively small region. 

The skewness of the velocity divergence is an interesting quantity to characterize the 
non-Gaussianity induced by gravity and is expected to constrain the density parameter with small 
theoretical ambiguities as long as primordial fluctuations are Gaussian distributed. But according 
to the present analysis we cannot determine the skewness of this field with an error less than 30% 
if our survey depth is not as deep as ~ 200/i _1 Mpc. 

In the previous Letter (Seto &; Yokoyama 1998) we have shown that the peculiar velocity 
field suffers from much larger sample variance than the linear density field because the former 
depends on small wavenumber modes much more strongly. On the other hand, the peculiar 
velocity divergence field discussed here has the same spectral dependence as the density field 
in linear theory (eq. O). Hence the large relative sample variance we have encountered in the 
present Letter is entirely due to the fact that it is nonvanishing even in linear theory whereas the 
expectation values of the higher-order cumulants become nonvanishing only after nonlinear effects 
are taken into account. 

This work was partially supported by the Japanese Grant in Aid for Science Research Fund 
of the Monbusho Nos. 3161 (N. S.) and 09740334 (J. Y.). 
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Fig. 1. — Relative sample variance E(X, V)/| (X) \ as a function of the patch radius L with 
different density parameter Qq. Thick lines represent the cumulants (solid: second-order, short- 
dashed: third-order, long-dashed: forth-order) of the velocity divergence field 9{x) smoothed with 
a 12/i _1 Mpc Gaussian filter. Thin lines represent the counterparts for the density contrast field 
5(x). Two lines for the second-order moments of 9{x) and 5{x) are identical. 
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